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ABSTRACT 

In this paper we present the derivation as well as the numerical results for all electro- 
magnetic form factors of the nucleon within the semibosonized Nambu-Jona-Lasinio (chiral 
quark soliton) model. Other observables, namely the nucleon mean squared radii, the mag- 
netic moments and the nucleon-A splitting are also computed. The calculation has been 
done taking into account the quark sea polarization effects. The final results, including ro- 
tational 1/N C corrections, are compared with the existent experimental data and they are 
found to be in a good agreement for the constituent quark mass of 400-420 MeV. 
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1. Introduction 



In the last years the most challenging problem in elementary particle physics seems to 
be solution of QCD in the low energy region. The main difficulties are due to the non- 
perturbative effects caused by growing effective coupling constant of the fundamental theory 
in the low energy limit. This prevents one from using the well-known main tool of theoretical 
physics — the perturbation theory. Because of this, the most intriguing features of QCD - 
confinement and chiral symmetry breaking - still remain conceptual and practical problems. 

The above mentioned obstacles have initiated an increasing interest among the physi- 
cists in non-perturbative methods and effective low-energy models of hadrons. The effective 
models are expected to mimic the behaviour of QCD at energies below ~ 1 GeV (confine- 
ment and/or chiral symmetry breaking) and to reproduce experimental data in this region. 
In principle, those models could be related to QCD by integrating out gluonic fields and 
reparameterizing fermionic degrees of freedom (the bosonization procedure). In effect, the 
physically observed particles should appear as fundamental fields of the model. 

One of the first models of this type that has become very popular was the Skyrme model 
[I]. It is an example of a fully bosonized effective model. During many years of intensive 
investigations the Skyrme model has proven to be theoretically quite interesting, though 
comparison with experiment has not always been good. 

The simplest purely fermionic Lorentz invariant model with spontaneous chiral symmetry 
breaking is the Nambu-Jona-Lasinio (NJL) model [2] that has followed the example of the 
BCS model of superconductivity. It contains chirally invariant local four-fermion interaction 
terms and it is non-renormalizable in 4-dimensional space-time. The NJL Lagrangian in its 
simplest SU (2) form has the following structure: 



where G is the coupling constant, mo is the current quark mass and f are the Pauli matrices 
in the isospin space. Usually up and down quarks are assumed as degenerated in mass. 

The NJL model is generally solved after applying the well known bosonization procedure 
following Eguchi [3] to arrive at the model expressed in terms of physically observed particles' 
fields cr(x), n(x): 



£njl = q - mo) q + 2 G [ (^) 2 + (w^q) 2 } , 



(1) 



£>NJL = q I il^d^-g(a + i7VT-f 5 ) ]q- -/i 2 (a 2 + tt 2 ) + (m fJ, 2 /g)a . 



(2) 
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The new coupling constant fj 2 is related to the initial G by: G — g 2 /fi 2 - Here, g is the 
physical pion-quark coupling constant implying that n are the physical pion fields. The 
meson fields are constrained to the chiral circle to reproduce the pion decay constant: 

a 2 (x) + f 2 (x) = f 2 , (3) 

where f % = 93MeV is the pion decay constant. 

In fact, an equivalent effective chiral quark meson theory [4,5] can be derived from the 
instanton model of the QCD vacuum. 

In the chiral quark soliton model [5-7] based on the lagrangean (2) (frequently reffered 
simply as NJL model) the baryons appear as a bound state of N c valence quarks coupled to 
the polarized Dirac sea (qq). Operationally the baryon sector of the model is solved in two 
steps. In the first step, motivated by the large N c limit, a static localized solution (soliton) is 
found. It is done by solving the equations of motion in an iterative selfconsistent procedure 
assuming that the a and n fields have hedgehog structure. However, this hedgehog soliton 
does not preserve the spin and isospin. In order to describe the nucleon properties one 
needs nucleon states with the proper spin and isospin numbers. To this end the classical 
solution is semiclassically quantized making use of the rotational zero modes. In fact, such 
a cranking procedure was elaborated in [5]. Recently, within this scheme including Dirac 
sea polarization effects quite successful calculations for the nucleon moment of inertia [13], 
magnetic moments [14] and nucleon electric form factors [15] have been reported. 

There are several reasons that the model described by (5) is considered as one of the 
most promising effective theories describing low energy QCD phenomena. First, the model is 
the simplest quark model which provides mechanism for spontaneous breaking of the chiral 
symmetry — the very basic feature of QCD. The mesons appear as excited states of quark- 
antiquark pairs. The philosophy behind this approach is based on the hypothesis that the 
chiral symmetry breaking and confinement are relatively independent mechanisms and they 
can be investigated separately [5]. 

Second, the self-consistent solitonic solution (hedgehog) has been found to exist in the 
large N c limit of the model within the physically acceptable range of parameters [7]. The 
chiral soliton provides a good description of the nucleon and gives us possibility to take 
into account vacuum polarization effects from the Dirac quark sea. Calculations done in 
the semiclassical quantization procedure yield quite reasonable results for the vacuum and 
nucleon sectors [8-10]. In particular, quantities like the nucleon radii, E-terms, axial vector 
coupling constant, A-nucleon splitting, splittings within the octet and decuplet are well 
reproduced. 
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Third, there are various hints that the NJL-type Lagrangian can be obtained [4,5,11,12] 
as a long wavelength limit of QCD. It should be stressed that the large N c limit plays a 
prominent role in those considerations. 

The aim of this paper is to compute the electromagnetic nucleon form factors within 
the SU{2) chiral quark soliton model based on the semibosonized NJL-type lagrangean 
with the vacuum polarization effects taken into account. The calculations will include the 
rotational 1/N C corrections, which really have been shown to be important for the isovector 
magnetic moment. Calculations of this sort give us possibility to calculate all electromagnetic 
properties of nucleons, as well as such quantities like the electric mean squared radii, the 
magnetic moments, the nucleon-A energy splitting, and the electric and magnetic charge 
distributions. 

The paper is organized as follows: we begin in Sec. 2 with formulation of the regularized 
model, introduction of the rotational zero modes and we derive expression for the electro- 
magnetic current treating the angular velocity as a perturbation. In Section 3 we derive the 
analytical expressions for the form factors including terms up to the linear order in angular 
velocity. Section 4 is devoted to presentation of numerical results. Summary and discussion 
is given in the last section. 



2. The regularized action and current 

In Minkowski space the generating functional of the model written in terms of the quark 
and meson fields, a(x),ir(x) has the form: 



Z'njl = J 'DqDqDv'D'K exp i J d 4 xC' NJL (x) 



(4) 



The sources are not explicitly included. 

Since the NJL model is a non-renormalizable theory, a regularization scheme of an ap- 
propriate cut-off A is needed to make the theory finite. To this end we will work in Euclidean 
space-time (the conventions used are given in the Appendix). 

After Wick rotation (transformation to the Euclidean space-time) integrating out the 
quarks in (4) one gets the effective action split in fermionic and mesonic parts: 

S eff = S eff + S ^ff + S eh , (5) 
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The fermionic part of the action includes the quark loop contribution in the presence of the 
external meson fields (a, n): 

S? ff = -8j>H0E)- (6) 
The Dirac operator Ip E (in Euclidean space-time) has the form: 

JP E = 74 [d A - h(U)} . (7) 

with the one particle Hamiltonian operator h E defined as 

h E (U) = i [-i lA lkdk + 74 MU\ . (8) 

The constituent quark mass M = gf^ and 

U = ^-(a + $757*7?) (9) 

represents the mesonic fields. It can be checked by inspection that the Hamiltonian is 
Hermitian: {fiEy = h E - 
The mesonic part 

S ^ff = Y jd A x E {a 2 + ^) ) (10a) 

and 

S b J ff = J d^x E a , (106) 

includes an explicit symmetry breaking term (10b). 

In order to compute the electromagnetic current matrix elements we include the electro- 
magnetic coupling in a minimal way and we take into account the rotational zero modes. 

We start with the Dirac operator^ with electromagnetic potential 

lp E = li [d 4 - h(U) - iA 4 Q + i l4lk A k Qj , (11) 

The quark charge matrix Q = ^1 + ^r 3 . 

In the next step, we introduce the rotating soliton: 

U(x, x 4 ) = R{x 4 ) U(x) R\x 4 ), R e SU 2 , R ] R = 1. (12) 

where the meson field U(x) corresponds to the stationary point (minimum) of the effective 
action (5). In the Euclidean space-time the rotation matrix R(x 4 ), can be represented as: 

R(x 4 ) = exp(-ix 4 £1) , ft = \r A Q A and Q E = -i R) R , (13) 
where as usually the summation convention over repeated indices is assumed. 
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In body-fixed frame the Dirac operator of the rotating soMton lp E (U) can be expressed 
in the form: 



where Ip E is defined by: 



ip E (U) = 74 \d A - h(U) iA A R ] QR + mikA k R ] QR 

The operator Q is Hermitian in both, Euclidean and Minkowski space-time: 



(14) 



(15) 



n E = n E , ti M = n M . (16) 

Taking into account (A.l) and (16) we can write down the Hermitian conjugate to (15): 



-d 4 - h{U) -iVL + iA\ R ] QR + iiAlkA\ R ] QR 



74 



;i7) 



where the collective coordinate Qm is will be quantized according to the canonical quanti- 
zation rule, which in Minkowski space-time are: 



A 



I 21 



(18) 



Here T A are the SU (2) group generators. The moment of inertia / has been computed 
numerically in [13]. 

Now, using the ansatz (14) the generating functional (4) can be re-expressed in terms of 
the rotational matrix R(t) in Euclidean space-time as 



Z'njl[\] = J VR(t) exp[-S eff [A^R(t),ny 



(19) 



where Q is a function of R = R(t) (see eq. (13)). Having the explicit form of the Dirac 
operator (15) we can go back to the action (5). Since we treat the meson fields classically 
(i.e. at the 0-loop level) the only non-trivial part of the action is the fermionic part (6). 
By construction (see [5,7]) it is splitted in a natural way into valence and Dirac sea parts. 
The valence part is finite and it will not be regularized. In this case the calculation is 
straightforward and the detailed derivation can be found e.g in [5]. Hence, from now on we 
will concentrate on the sea contribution and we will restrict our discussion to that part of 
the full action. 

As we can see from (15), (17), the operator (15) is non-Hermitian. Hence, in general, the 
action (6) will have the real and imaginary part. To make clear distinction we write down 
both parts explicitly: 
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S? ff = ReS? ff + ImS? ff , (20) 
ReSL f = -\s V H$$] , (20a) 



ImS[ ff = -±8p\n[9>/$] , (206) 

The imaginary (anomalous) part is finite and does not need regularization. In addition, 
any regularization of this part of the action would break several important features of the 
model related to the anomalous structure. On the other hand, the real part is infinite and 
a regularization is necessary. 

In order to compute the nucleon form factors we have to evaluate the expectation value 
of the electromagnetic current j^ m . To this end we take into account the standard definition 
given by: 



{j em m ^ 1 SZ' NJL [A,} 



Z SA^(0) 



i r 5S eff 



a„=o 



exp[-S e// ], (21) 



A,=o Zo J ~ ' SA^O) 
where the generating functional Z' NJL is given by (19). 

As it was mentioned above the real part of the action diverges. To calculate the expression 
(21), we introduce the Schwinger proper-time regularization of the action. In general, for an 
operator A its proper-time regularized form reads [17]: 



hiA = - lim / — e~ TA . (22) 
A ^°° Jl/A 2 T 

The Schwinger method is especially convenient for the logarithmic expressions resulting from 
the fermionic determinant. Applying (22) to the action (6) we obtain: 

S eff = yE d 4 x E Tr^Tr T I —tf n (x E ) exp[-r^^] ^ Xe ) , ( 23 ) 

where ij) n (x) is any complete set of eigenf unctions. Actually, in our calculations we shell use 
for ip n (x) eigenfunctions of the Hamiltonian h. For the meson fields a, n we substitute the 
static selfconsistent hedgehog solution U (x) of the chiral soliton [7] form: 



a(x) = cr(r), tt(x) = rir(r) , (24) 

restricted to the chiral circle (3). Numerically, we will follow the method of Ripka and 
Kahana [16] and will use the eigenfunctions of the Hamiltonian h to compute all matrix 
elements (see also Section 4). 
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Now, our task is to evaluate the current (21) with the regularized action (23). The diffi- 
culty lies in the fact that the expression (23) contains rather complicated operator exponent 
in the integrand. However, because afterwards we take derivative over the fields and put 

= 0, the only nonzero contribution will come from the part of the integrand linear in A^. 
We shall treat the cranking parameter f] ~ as a small perturbation that is consistent 
with the large N c limit philosophy behind the NJL model. Hence, we can neglect terms Vt 2 
and higher*. Also, we should be careful in performing the VR integration, as in general 
the operators Q and R does not commute and one must use the proper time ordering (see 
Sect. 4 and [18, 19]). In the next step we expand the integrand in (23) in terms of A^ and 
fl To this end we apply the Feynman-Schwinger- Dyson (FSD) expansion for the operator 
exponent: 

e A+B = e A + f 1 dae aA Be (l-a)A + ^ ^ P da £ aA g e (3A g e (l-a-0)A + (35) 

Jo Jo Jo 

In our case, for the operator A we substitute the part of our operator exponent (—tW), 
where: 



W=ip E ^Ip E , (26) 

that does not contain A^ and f2. This part of W we will call Wo- On the other hand, 
for the operator B we take the remaining terms of the operator W and we will call them 
Wi, W2, W5. The first type of terms (Wi) is defined as the sum of all terms in the exponent 
that are linear in A^, W2 ~ A^ x f2, W3 ~ f2, W4 ~ f2 2 and W§ ~ A^. The W§ contribution 
to (21) is exactly zero, while the W4 contribution we neglect as being small (quadratic in 
Q). Hence, the operator B will include the terms W1-W3 only. Finally, this gives us the 
following expressions for the terms Wq-Ws of the operator W (26): 



W = + cu 2 + h 2 , (27a) 
Wi=- (A 4 + 4) R ] QRuj + (A k - A\) 7 47fc R ] QRuj + 
+ i (h A 4 R^QR - A\ R^QR h) - 

- (h 141k A k R^QR + 7 47fc A k R^QR h) , (27b 1 ) 



* Actually, the terms ~ f2 2 must be taken into account to evaluate the moment of inertia 
(see [5, 13] for details). 
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W 2 = - (A 4 Q R ] QR + A\ R^QR SI) + 747^ (A k n R^QR - a\ R^QR Q) + 

+ A 4 [Q, R^QR] - l4lk A k [ft, R^QR] , (27c') 
W 3 = + 2VLu -i[h,tt] . (27d) 

Taking into account that the electromagnetic potential is Hermitian we get for the W\ and 
Wi the simpler form: 

Wi = + % [h, A 4 R)QR] - i {h, 747fc A k R^QR} - 

-2A 4 R^QRlo, (27b) 
W 2 = -2A 4 VtR^QR . (27c) 

where [ , ] and { , } denote commutators and anticommutators, respectively. Also, the sum 
of all perturbative terms that will contribute the electromagnetic current (21) we will denote 
by V: 

V = W 1 + W 2 + W 3 . (28) 

In addition, the part of V containing (linear) u;-terms will be called uj V\ and the u;-indepen- 
dent part will be denoted by Vo- Hence, the full perturbative term reads: 

V = Vq J ] + Vi u 1 . (28 ; ) 

To obtain (27) we have computed the trace over the Euclidean time in the cu-space, the 
Fourier conjugate to x 4 . To this end we have substituted in (23): 



/ 



f /.n /v.', j I ^-fi+icj). (29) 



Also, the formulae (A.4-A.5) from the Appendix have been used. 

Applying (27) and (25) we should have in mind that there are, in general, two types 
of terms in (25): linear and bilinear in B. Hence, the products of the terms (27b-27d) can 
contribute as well. Also, those products should be linear in f2 and A^. This implies that the 
only rV-term contributing to the second order FSD expansion is of the form: W\ x W3 - 
the first one linear in the potential A^ and the second one linear in the cranking parameter 
Q. As we will see later on, those terms will contribute to the electric isovector and magnetic 
isoscalar parts of the form factors. 
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As the above considerations were limited to the regularized part of the action a comment 
concerning the non-regularized parts is in order. In general, there are two types of non- 
regularized contributions: the valence parts (for all form factors) and the sea parts of both 
isoscalar form factors. To compute these quantities we can start directly from the non- 
regularized action (6) as it was done in [5]. The electromagnetic current reads: 



/d 1 
2^ ^ n{x) lo + ih - n R]QR 74> Mx) ' (30) 
n* - it 

where the integral over oo stands for the trace in Euclidean time according to (29). Further 
calculations can be easier done, without using the expansion (25), just expanding (30) in 
terms of (small) Vt. The same results can be obtained by imposing the A 2 — > oo limit on our 
regularized expressions (using formulae (A. 2) of Appendix). 

The second remark concerns our practical calculation of the imaginary part of the action. 
One can notice that taking reversed Hermicity (anti-Hermitian) electromagnetic potential 

in (27b', 27c') the imaginary terms are moved to the real part. With this trick and 
the A 2 — > oo limit we can calculate all form factors from (27) within our scheme and not 
referring to the nonregularized action (6). 

Now, taking into account (21), (25) and (27-28) we can write down the operator exponent 
contributing to the electromagnetic current as: 



M x ) = 



Seff[A,} 



5A^ 



/ — exp[-rW] ~ -t / — \ e ™ (n\V\n) + 
J-oo 2tt J_ 00 2tt ^ 

+ r 2 / + °°^y e -™ 2 \(n\V\m}\ 2 ['dp C~ P da e -M)4 e -r^ m ? 

J-oo 2tt Jo Jo 



m,n 

where the terms quadratic and higher in Q are neglected and because of the derivative in 
(21) only the terms linear in potentials are present. To compute the functional (space-time) 
trace we use the Hamiltonian (8), like for the action (23). Here, the vectors \n) denotes the 
eigenfunctions i/) n (x) corresponding to the eigenvalues e n . 

We cannot perform the tu-integrals directly as the operator V is u;-dependent. However, 
we can simplify (31) taking into account that the term |V^, n | 2 is symmetric in indices m, n 
and performing trivial a- and /5-integrals. After that the integration over uj is performed 
and we get for the electromagnetic current the following expression: 



10 



Jem 



N c f°° dr_ f 



4Vtt Ji/a 



-re 2 / |^0 



n + 



_ 2 _ 2 



m,n 



e 2 -e 2 



5^0 



(n\V [A^ = 0]\m){m\—f-\n) + 



5A, 



(32) 



+ — (n\V 1 [A fl = 0}\m)(m 
It 



5Vi 
SA U 



\n) 



where the notation V\ [A^ = 0] means the sum of all terms linear in u which do not contain 
Afj,. The formula (32) together with (27) will be used in the next Section to obtain the 
electromagnetic form factors' formulae. 



3. Computing the form factors 

In this point, applying the FSD expansion (25) we are ready to compute the electromag- 
netic form factors directly from the regularized action (23) which includes the fermion loop 
corrections. The nucleon electromagnetic Sachs form factors are defined by: 



(N(p)\f m (0)\N(p')) = G E (q 2 ) , (33a) 
(N(p)\jf m (0)\N(p')) = ^ G M (q 2 ) i e ikl q k (N\a l \N) , (336) 

where \N(p')) is the nucleon state with proper spin and isospin quantum numbers. The 
electromagnetic current j^ m (a;) is in the Minkowski space-time and M.^ is the nucleon mass. 
It is clear from (33) that to compute the electromagnetic form factors we need the current 
evaluated in the previous Section. The isoscalar and isovector parts of the form factors are 
defined by: 

G e ,m = 2 g e\m + ^3 Gem ■ ( 34 ) 

The Sachs form factors (33) are appropriate for the non-relativistic limit (q 2 <C M-n) that 
is consistent with the large N c approximation. 

The isoscalar form factors (electric and magnetic) come from the scalar part of the quark 
current matrix Q (~ 1) while the isovector form factors is from the triplet part (~ r 3 ). The 
terms linear in Q contribute to the electric isovector and to both magnetic form factors. As 
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can be easily checked directly from (27) and (32), the contribution from the real part of the 
action to both the isoscalar electric form factor and the isoscalar magnetic one is exactly 
zero. This means that these form factors originate from the imaginary part of the action 
and they will not be regularized. Table 1 lists terms contributing to a given form factor. In 
the last column we explicitly indicate which form factor should be regularized. 

Table I. Terms of the effective action contributing to different electromagnetic form factors. 



quantity 


\ 




Q 


2°-order FSD 




regularized 


r<T=i 
r<T=i 

M 


A± 
A± 
A k 
A k 


n° 
n 1 


~ i 
~ i 

~T 3 


NO 
YES 
YES 
NO + YES 


NO 
YES 
NO 
NO 


NO 
YES 
NO 
YES + NO 



The electric scalar form factor can be immediately extracted from (33a) and (30). Since 
only the isoscalar part of the operator Q contributes the rotational matrices R(t) vanish 
from the action. Hence, the D1Z integration in (19) becomes trivial. Also, in this case the 
contribution from the terms linear in Q is exactly equal to zero. Finally we arrive at: 

G£=V) = / d 3 x e-& ^ { £ limnix) ~ 

J ^ neVAL 

(35) 

neALL ' 

where val denotes the valence whereas all stands for all eigenstates of the Hamiltonian 
including the valence one as well. The first term in eq.(35) comes from the valence quarks 
and the second term is due to the quark vacuum polarization. As has been stated previously, 
the same result can be obtained by substituting ant i- Hermit ian fields in (27b'), using 
(32), and taking the limit A 2 — > oo. 

As the next step, we compute the electric isovector form factor. In this case, the regu- 
larization is necessary and the second order terms in the FSD expansion have to be taken 
into account. Contribution to the first order come from the terms of Wi containing both A± 
and Q. The second order includes terms from W\ and W3. 
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In the first order FSD term we insert unity: J2 m il>m(x)' l Pm{ x ) an d we get identical overall 
isospin-space factor of the form: 

ji m ~ ^ n Q^ m ^ m Rr 3 R^ n , (36) 

to which the Fierz identity (A. 6) is applied. Then, we express Q in terms of the rotational 
matrix R using the identity: 

Q = —iRR} = -iR}R (37) 

and (A. 7). 

Finally, from the quantization rule for the collective coordinate Q (18), where the moment 
of inertia I has been computed numerically in [13], we arrive at the following expression for 
the electric isovector form factor: 



G T E =\f) 



d 3 x e~ iqx 



-N c 
61 



neVAL 
m£ALL 



+ 



1 v ^ r°° dr_ 



en + e r 



2 2 

+ T 



X 



(38) 



The formula (38) for q 2 = is equivalent to the quotient of the moment of inertia by 
itself and: G 1 ^ =l (q 2 = 0) = 1, as it should be. 

We proceed to compute the magnetic form factors. We start with the isoscalar one. In 
this case, the contribution comes from the imaginary part of the action and the regularization 
is not necessary. The quark current matrix Q contributes via the isoscalar term, ~ | 1, 
diagonal in the isospin space. Hence, the rotation matrices R(t) cancel each other and, like 
for the electric scalar form factor, there is no time ordering problem for the DR integration 
in (19). However, there is still a nonzero contribution from terms linear in Q. Starting 
from the reversed Hermicity of and the regularized (real) action term (27c') we find 
that contribution to the first order FSD expansion comes from the anticommutator term: 
747^ Aj, {Q, RQR)} = |747^AfcQ. It contributes to the electromagnetic current as follows: 



k 

Jem 



(1) 



120FJ i 

A z 



dr 



(n|747 fc ^|n) . 



(39) 
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The second order FSD contribution comes from the terms 3^747/^ — ^[/i, 747^^] (in 
W\) and 2VLuj — i[h, Q] (in W3). Taking into account (32) and (39) together with (13) and 
(18) we obtain the space components of the isoscalar electromagnetic current in Euclidean 
space-time in the following form: 



fk m = 



N r _ f°° rfr v 



77 



_ 2 _ 2 

c p T€ n _l_ (z p T€ m 



1 e" re * 

+ o 



X 



(40) 



x / A(iV|r A |iV)^t(f)747fc^m(x)^(y)r^ n (y) 



Now, we have to perform transformation to the Minkowski space-time. Also, one has 
once again to resolve (33b), this time with respect to the form factor G^f°. To this end we 
apply (A. 9) and the result for the sea part of the magnetic isoscalar form factor reads: 



G T=0,„2 



r M 



-I 



d s x e~ iqx 



x 



-N C M N f 


d 3 ye ijA % x 


f°° dr 




2 2 
< fc me 






e n "1" e m 



+ 



le" 



e 2 -e 2 



X 



(41) 



X ^l(xholi^m(x) i>]n{y)T A i> n (y) . 



t 



To obtain the full form factor we must add to (41) the valence contribution and we must 
take the non-regularized limit: A 2 — > 00. For the last step we use formulae (A. 2) to find 
that: 



lim 

A 2 ^oo 



1/A 2 Vt 



_ 2 _ 2 

<r p T6 n _|_ , p T€ m 



+ 



(42) 



7T 



[ sign(e n ) - sign(e r 



Finally, we get for the magnetic isoscalar form factor the following expression: 
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n^VAL tn ~ €m 

meALL 



(43) 



1 . f v i)n{x)lQlii)m{x) 1ptn{y)T A 1p n {y) 

+ 2 2.^1 sl S n l e nJ 

rieVAL 
mGALL 



where, as usually, summation over repeated indices i,j,A is assumed. The first part gives 
the valence and the second part the sea contribution to the form factor. 

For the magnetic isovector form factor the sea contribution comes from the isotriplet 
part of the quark current matrix Q. Because of this we cannot get rid of the rotation matrix 
R(t) in the action. 

First, let us compute contribution from the first order FSD expansion (0-order in Q). 
This term is the anticommutator — 747^ Afci?r 3 in (27b). This implies for the elec- 
tromagnetic current: 



= E If ^ €ne ~ Te "H747 fc r 3 |n> , (44a) 



n " A 2 



where we are still in the Euclidean space-time. Taking into account (A. 8) and 70 instead of 
74 we get: 

_ at roo 1 

fk m = 7tt Tr (R^R^r A )J2 ~r^n e—" <n| 7 o7*r» , (446) 

For the term (44b) regularization is necessary, as it comes from the real part of the action 
and is infinite. 

For the linear order in Q we have to take into account non-commutativity of the isospin 
operator in the collective variable (18) with the finite rotation (Wigner) matrix. This leads 
to the extra term, linear in the cranking parameter The observation has been done in [18] 
in the context of the axial coupling constant g^, and the proper treatment of such effect for 
the qa and the isovector magnetic moment has been given in [19]. 

Following the method of [19] we get the additional, linear in f2 correction to the isovector 
electromagnetic current in the form: 
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jr A (n 1 ) = N e [n c ,\Tr(BtT A RT B )] l d "y 

™r\r at J 



3 ^U^holkT B iJ m (x)ijln(y)T C ij n (y) 



m£VAL 
neALL 



(45) 

Here, a remark concerning regularization of the term (45) is in order. This expression is nu- 
merically finite and, as has been shown in [19], it is finite in the gradient expansion. Hence, 
the regularization is not necessary. Also, one should have in mind that the main contribu- 
tion to the isovector magnetic current comes from the valence part, the sea contribution is 
relatively small and the formal regularization of that contribution gives a small difference 
with our result. 

To calculate the form factor we have to resolve the equation (33b) with respect to G^f 1 
using the spherical basis in which the rotation matrices R(t) are related to the Wigner func- 
tions. Our current (44b) has the general index structure of the form: = X A tkAili F(q). 
Comparing this structure with eqs. (33b, 34), using (A. 9) we obtain the final expression for 
the magnetic isovector form factor in the following form: 

G T M =\ q 2 ) = [ d\ e-& e ijA % \ £ w ^ n(f) _ 

J 6 q ^neVAL 

- ttt fi ~T €n e ~ T£ " ^n( f )707i^Vn(x) + , , 

1 f J3 ri(£)707* T B ll>m(x)A(y)T C ll> n ( 



+ h tABC S J d 



n>val 
m<val 



The proton and neutron form factors are defined, respectively, as sum and difference of 
the isoscalar and isovector form factors: 



2 

1 
2 



Gem ~ o [ Ge,m ~ G e ,m ] > 



and the numerical results for them are presented in Section 4. 



4. Numerical results 
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To compute observables we use the finite quasi-discrete basis and numerical method of 
Ripka and Kahana [16] for solving the eigenvalue problem. The Hamiltonian h is taken in 
a spherical box of a large radius D. The basis is made discrete by imposing a boundary 
condition at r — D. Also, it is made finite by restricting momenta of the basis states to be 
smaller than the numerical cut-off K max . Both quantities have no physical meaning and the 
results should not depend on them. The typical values used are D ~ 20/M and K max ~ 7M. 

In addition, all checks concerning the numerical stability of the solution with respect to 
varying box size and choice of the numerical cut-off have been done and the actual calculation 
is completely under control. 

The actual calculations are done by fixing the parameters in meson sector in the well 
known way [7] to have /„• = 93 MeV and = 139.6 MeV. This leaves the constituent 
quark mass M as the only free parameter. 

The proton and neutron electric and magnetic form factors are displayed in figs. 1-4. The 
theoretical curves resulting from the model are given for the following four values the con- 
stituent quark mass: 370,400,420 and 450 MeV. The magnetic form factors are normalized 
to the experimental values of the corresponding magnetic moments at q 2 = 0. With one 
exception of the neutron electric form factor (fig. 2), all other form factors fit to the experi- 
mental data quite well. The best fit is for the constituent quark mass around 400-420 MeV. 

As can be seen the only form factor which deviates from the experimental data is the 
neutron electric form factor and this requires some explanation. Obviously, this form factor 
is the most sensitive for numerical errors. According to the formula (47b) the form factor has 
been calculated as a difference of the electric isoscalar and electric isovector form factors that 
were the direct output of our code. Both form factors were of order 1.00 and calculated by the 
code with rather high accuracy of about 1%. However, the resulting neutron form factor has 
the proper (experimental) value of order 0.04, i.e. about 4% of the value of its components. 
This means that a small numerical error for one of the components can be enhanced by 
a factor 50. Hence, small numerical errors together with the applied approximations (like 
the 1/N C approximation behind the model and neglecting the boson-loop effects as well as 
higher order fermion loops) are strongly magnified resulting in a considerable deviation from 
the experimental data for momentum transfers above 100 MeV j c. 

As the next step, one can compute other electromagnetic observables: the mean squared 
radii, the magnetic moments and the nucleon-A splitting. The static nucleon properties, in 
particular the charge radii and the magnetic moments can be obtained from the form factors: 
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^ = V) 2 • (49) 

For the quark masses 370,420 and 450 MeV and pion mass m^ = 140MeV the calculated 
values are presented in table 2. For comparison, in table 3 we give the theoretical values for 
the same quantities but with the physical pion mass set to zero. 

Table 2. Nucleon observables computed with the physical pion mass. 



Quantity 


Con; 

370 MeV 


stituent Quark 1 

420 MeV 


Vlass 

450 MeV 


Exper. 


total 


sea 


total 


sea 


total 


sea 


< r 2 > T =o [fm 


21 


0.63 


0.05 


0.52 


0.07 


0.48 


0.09 


0.62 


< r 2 > T =l [fm 


21 


1.07 


0.33 


0.89 


0.41 


0.84 


0.45 


0.86 


< r 2 > p [fm 2 




0.85 


0.19 


0.70 


0.24 


0.66 


0.27 


0.74 


< r 2 > n [fm 2 ] 


-0.22 


-0.14 


-0.18 


-0.17 


-0.18 


-0.18 


-0.12 


Ht=o [n.m] 


0.68 


0.09 


0.62 


0.03 


0.59 


0.05 


0.88 




n.m. 




3.74 


0.96 


3.53 


1.07 


3.47 


1.17 


4.71 


Hp 


n.m. 




2.21 


0.53 


2.08 


0.55 


2.03 


0.61 


2.79 




n.m. 




-1.53 


-0.44 


-1.46 


-0.52 


-1.44 


-0.56 


-1.91 


[Vp/Hnl 


1.44 




1.42 




1.41 




1.46 


M A - M N [MeV] 


213 




280 




314 




294 


9A 


1.32 


0.15 


1.29 


0.20 


1.29 


0.23 


1.26 



The chiral limit m^ — > mostly influences the isovector charge radius. In fact, as 
illustrated in fig. 5, where the isovector charge radius is plotted vs. the box size D, it 
grows linearly with D and diverges as D — > oo. Because of this that quantity (and the 
derivative quantities) are not included in table 3. The other quantity strongly influenced 
by the chiral limit is the neutron electric form factor. For the m n — > the discrepancy 
from the experiment is larger by almost a factor two than in the case m n ^ 0. The other 
quantities differ in the chiral limit by about 30%. The comparison of values in the two tables 
indicate that taking the physical pion mass gives us a best fit with a much better agreement 
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with the experimental data. To be particular, in the calculations with zero pion mass the 
electric observables suggest a high value (M ~ 450 MeV) while the magnetic ones indicate 
M ~ 370 MeV which is not the case with physical pion mass. In addition, we observe much 
larger contribution from the sea effects, about 50% of the total value. 

Table 3. Nucleon observables computed with the zero pion mass. 



Quantity 


Co 

370 MeV 


Qstituent Quark 1VI 

420 MeV 


ass 

450 MeV 


JJAL/Cl • 


total 


sea 


total 


sea 


total 


sea 


< r 2 > r =o [f™ 2 




0.88 


0.20 


0.66 


0.26 


0.61 


0.23 


0.62 


(J,T=o [n.m. 




0.66 


0.07 


0.59 


0.09 


0.57 


0.09 


0.88 [ 




n.m. 




4.84 


1.98 


4.50 


2.09 


4.29 


2.16 


4.71 


fi p 


n.m. 




2.75 


1.03 


2.54 


1.09 


2.43 


1.13 


2.79 


fi n [n.m. 




-2.09 


-0.96 


-1.95 


-1.00 


-1.86 


-1.03 


-1.91 






1.31 




1.30 




1.31 




1.46 


M A - M N [MeV] 


221 




261 




301 




294 


9A [MeV] 


1.34 


0.24 


1.29 


0.15 


1.25 


0.00 


1.26 



The results of table 2 (m n = 140 MeV) again indicate the value ~ 400-420 MeV for 
the constituent quark mass, in agreement with the conclusion drawn from the form factor 
curves. The same has been suggested earlier [13-15], where smaller number of observables 
has been taken into account. With the exception of the neutron electric squared radii, to 
which remarks similar to the case of the neutron electric form factor are applicable, the 
contribution of the valence quarks is dominant. However, the contribution of the Dirac sea 
is non- negligible, it lies within the range 15 - 40%. 

One can notice, that the numerical results for the nucleon-A mass splitting (Ma — Mn), 
the mean squared proton, isoscalar and isovector electric radii and the axial coupling constant 
(<m), as well as the proton electric and magnetic and neutron magnetic form factor, differ 
from the experimental data by no more then about ±5%. Finally, for the magnetic moments 
we have got results 20-25% below their experimental values. Despite of this underestimation 
of both magnetic moments we have found surprisingly good result for the ratio \x v j \i n which 
is far better than in other models. The agreement with the experimental value is better than 
3% for the constituent quark mass 420 MeV . So good result has been obtained neither in 
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the Skyrme, nor in the linear chiral model. 

The isoscalar and isovector electric mean squared radii are shown in figs. 6-7 as functions 
of the constituent quark mass. The same plot but for the experimentally measured quantities: 
the proton and neutron electric charge radii, is given in fig. 8. In these plots the valence and 
sea contributions are explicitly given (dashed and dash-dotted lines, respectively). As could 
be expected, for the isoscalar electric charge radius the valence part is dominant (about 
85%), due to the fact that there is no contribution from the terms linear in the collective 
coordinate Q. This is not true for the isovector electric charge radius, where the sea part 
contributes to about 45% of the total value (fig. 7). Also, this effect can be seen from the 
proton and neutron charge radii (fig. 8) which are linear combinations of the isoscalar and 
isovector ones. For proton the sea contribution is about 30%. However, for the neutron 
charge radius the sea part is dominating and the valence contribution is negligible. 

In addition, in fig. 9 we plot the magnetic moment density distribution for proton and 
neutron for the constituent quark mass M = 420 MeV. The sea contribution becomes non- 
negligible for distances greater than 0.5 fm. Due to the relatively large tail contribution 
to the magnetic moments the sea contribution to these quantities is about 30%. If one is 
interested in the proton and neutron charge distribution the reader is referred to ref. [15]. 

Both the magnetic moments and the electric radii are rather sensitive to the tail be- 
haviour of the corresponding densities. In actual calculations of these quantities the densi- 
ties are integrated with extra factors r for the magnetic moments and r 2 for the radii that 
enhances the role of the tail. This leads to the strong enhancement of numerical errors and 
partially it explains the deviations. 

5. Summary and conclusions 

Our numerical results support the view that chiral quark soliton model of the Nambu- 
Jona-Lasinio type offers relatively simple but quite successful description of some low-energy 
QCD phenomena and, in particular, of the electromagnetic properties of nucleons. For the 
first time the magnetic form factors are calculated. We have obtained quite good results for 
the electromagnetic form factors, the mean squared radii, the magnetic moments and the 
nucleon-A splitting reported in figs. 1-4 and tables 2, 3. 

Using = 93 MeV and = 139.6 MeV and a constituent quark mass of M — 
420 MeV the isoscalar and isovector electric radii are reproduced within 15%. The magnetic 
moments are underestimated by about 25%, however, their ratio fj, p / /i n is almost perfectly 
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reproduced. The g-dependence of G p E (q 2 ), G p M (q 2 ) and G 7 l I (q 2 ) is very well reproduced. 
The neutron form factor G E (q 2 ) is by a factor of two too large for q 2 > 100 MeV 2 . One 
should note, however, that G E (q 2 ) is more than an order of magnitude smaller than G p E (q 2 ) 
and as such it is extremely sensitive to both the model approximations and the numerical 
techniques used. Altogether it is fair to say that the NJL in SU (2) reproduces reasonably 
well the electromagnetic properties of the nucleon. Here, it turns out that the agreement is 
noticeably worse if the chiral limit, m n — > 0, is used. This can easily be understood since the 
pion mass determines the asymptotic behaviour of the pion field. Altogether we conclude 
that the bosonized NJL model chiral quark soliton model based on a bosonized NJL type 
lagrangean is quite appropriate for the evaluation of nucleonic electromagnetic properties. 
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Appendix 



In Euclidean space we use the metric tensor's signature transformed from (H ) to 

( ) (i.e. the Euclidean metric tensor gjjj v = —5^) and the general formulae to perform 

the transformation from Minkowski to Euclidean space are: 

E , • M 4 -o E , M 

a 4 = +ia , a E = -ia M , a k = +a k , 

„,E | -M „,E | „,M 

74 = + ? 7o » 7fc = +7fc , 

(Al) 

{Tf,7f} = -2^, {7jf,T#} = 0, 

(7?)* = -7?, (7 5 £ ) t = +7 5 £ = 7 5 M • 

Here, the simple substitution rule for the limits of basic r-integrals are given to calculate 
the non-regularized formulae for observables: 

[°° dr 2 
lim / —j- a e Ta — > sign(a) , (A2a) 

/■°° r/r 2 

lim / _ e -™ _>-27r|a| , (A26) 

where a is a real number. In addition, we can express the sign(x) function in terms of the 
step function as: 

sign(:r) = G(x) - Q(-x) = 1 - 20(-x) . (A3) 

To obtain the formulae (27) we have applied simple identity for the time derivative of 
the operator RQR) ', implied by (13): 

d A (RQR ] ) = i [RQR\ Q] . (A A) 

Also, one should have in mind that our potentials A^(x) are time- independent and we can 
use identity: 

<9 4 A = d A A k = . (A5) 

For the second order FSD expansion it is useful to apply the Fierz identity. If the 
matrices A, B: A= \A A t A , B = \B a t a , then we get: 

Tr (Ar A ) Tr (Br A ) = 2 Tr (AB) . (A6') 
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This can be applied to our overall isospin factor in the electromagnetic current as follows: 



(A6) 



J2 / dn s dn^l{y)^ m {y)^ m {x)B}T Z R^ n {x) = 

m,n 

J2\mvt a ) Tt(r^t 3 Rt b ) I dn s dn y ^t(y)T A ij m (y)ijl(x)T B MS) = 

m,n 

J2^Tr(Qr A ) Tr(R^r 3 Rr A ) f dn s dn y ^(y)r B ^ m (y)^l(x)r B M^ = 

m,n 

J2l^R f r 3 R) ! dn s d^i{y^T A ^ m {^l{x)T A M^) ■ 

m,n 

Further simplification is done using the identity: 

Tr (QRt 3 R^) = -i Tr (RR^Rt 3 R^) = -% Tr (R?Rt 3 ) = ft 3 , (A7) 
where we have used (37) and the obvious identity: 

A = 1 Tr (Ar A ) r A . (A8) 
2 

To resolve the magnetic form factor from the right hand side of eq. (33b) we use the 
identity: 

1 - 1 * Qj ( A r\\ 

1 = 2 tmBl emB i ~~^r ■ t A -y) 

This can be applied to our expression for the current (44b). 
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Figure captions 

Fig. 1. The proton electric form factor for the momentum transfers below 1 GeV. 

Fig. 2. The neutron electric form factor for the momentum transfers below 1 GeV. 

Fig. 3. The proton magnetic form factor normalized to one at q 2 = for the momentum transfers 
below 1 GeV. The normalization factor can be extracted from table 2. 

Fig. 4. The neutron magnetic form factor normalized to one at q 2 = for the momentum 
transfers below 1 GeV. The normalization factor can be extracted from table 2. 

Fig. 5. The isovector charge radius as a function of the box size D for m n = and m n = 
139.6 MeV. 

Fig. 6. The isoscalar electric charge radius as a function of the constituent quark mass M. The 
valence and sea parts are marked by the dashed and dashed-dotted lines. 

Fig. 7. The isovector electric charge radius as a function of the constituent quark mass M. The 
valence and sea parts are marked by the dashed and dashed-dotted lines. 

Fig. 8. The electric charge radii of proton and neutron as functions of the constituent quark 
mass M. The valence and sea parts are marked by the dashed lines. 

Fig. 9. The magnetic moment density of proton and neuron for the constituent quark mass M = 
420 MeV. 
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